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Abstract 

For the case of a simple face-centered-cubic pliotonic crystal of homogeneous 
dielectric spheres, we examine to what extent single-sphere Mie resonance 
frequencies are related to band gaps and whether the width of a gap can be 
enlarged due to nearby resonances. Contrary to some suggestions, no spectac- 
ular effects may be expected. When the dielectric constant of the spheres 
is greater than the dielectric constant of the background medium, then for 
any filling fraction / there exists a critical £c above which the lowest lying Mie 
resonance frequency falls inside the lowest stop gap in the (111) crystal direc- 
tion, close to its midgap frequency. If < the correspondence between 
Mie resonances and both the (111) stop gap and a full gap does not follow such 
a regular pattern. If the Mie resonance frequency is close to a gap edge, one 
can observe a resonance-induced widening of a relative gap width by « 5%. 
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I. INTRODUCTION 



Photonic crystals, i.e., dielectrics with spatial periodicity, have triggered much interest 
recently One can picture a photonic crystal as a periodic arrangement of dielectric 

scatterers. For example, its dielectric constant £:(r) equals Eg for r inside the scatterer and 
Eb otherwise. Similarly to the case of an electron moving in a periodic potential, a photon 
traveling in a photonic crystal encounters a periodicaly changing dielectric constant. As a 
result, a gap can open in the electromagnetic wave spectrum, independent of its polarization 
and direction of propagation In a given frequency interval the density of states (DOS) 
can either be reduced down to zero (photonic band gap) or enhaced with respect to its 
vacuum value. Such a change in the DOS affects many physical quantities. The most 
transparent is the change in the spontaneous emission (SE) rate of embedded atoms and 
molecules. This can be demonstrated already at a relatively low refractive index contrast 
(fs 1.15) 1^. Suppression of the SE may have applications for semiconductor lasers, solar 
cells, hetero junction bipolar transistors, and thresholdless lasers |jl]]. On the other hand, 
the enhancement of the SE is a way to create new sources of light for ultra-fast optical 
communication systems [Q. Unlike conventional (electronic) crystals, photonic crystals are 
essentially man-made structures and their parameters can be changed at will. There is a 
common belief that in the near future such systems will allow us to perform many functions 
with light that ordinary crystals do with electrons. 

Thus far, the main emphasis in the study of photonic crystals has been on calculating 
the band structure Let / be the scatterer filling fraction, i.e., the volume of the 

scatterer(s) in the unit cell per unit cell volume. Once / is fixed, the spectrum is only a 
function of the dielectric contrast 6 = max{Es/Eb,Eb/Es)- In our paper we pose the question 
of to what extent the single sphere resonance frequencies are related to band gaps, whether 
a gap width can be enlarged due to nearby resonances, and what other effects, if any, single- 
scatterer resonances may have on the band structure and properties of a photonic crystal. 

In the following, we focus on the case of a simple face-centered-cubic (fee) photonic crystal 
of homogeneous spheres. There are at least two reasons to consider this case. Firstly, 
single-sphere resonances, known also as Mie resonances |TD[, are well understood and an 



analytic solution exists for them. In each angular-momentum channel characterized by the 
angular-momentum number I a single sphere has an infinite number of Mie's resonances. 
The properties of Mie's resonances are discussed in many monographs (see, for example, 
||1 1|| ) and we only emphasize the following ones. The sharpness of a Mie resonance decreases 
with increasing ar^, where is the sphere radius, a = Uy/Fb/c, u is the frequency, and c is 
the speed of light in vacuum. On the other hand, reflecting the centrifugal barrier increasing 
with / as /(/ + 1), the sharpness increases with /. The spacing between resonances in the 
higher frequency range is accounted for by the resonance condition avs = m7r/2, where m 
is an integer (if = 1). 

Secondly, fee structures of homogeneous spheres are among the most promising candi- 
dates to achieve a full photonic gap at optical and near-infrared frequencies. Indeed, in this 
frequency range one often uses collodial systems of microspheres which can self-assemble into 
three-dimensional fee crystals with excellent long-range periodicity ||T2| , |T3i . This long-range 
periodicity gives rise to the strong optical Bragg scattering, clearly visible by the naked eye, 
and already described in 1963 |T^. Both the case of "dense" spheres {Eg > Eb) Il3| and "air" 
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spheres (e^ < Eb) [0] can be realized experimentally. 

The outline of our paper is as follows. In the next section, main features of the method 
that is used to calculate the spectrum of electromagnetic waves in periodic structures are 
discussed. We also discuss the differences between and similarities of the electronic and 
electromagnetic bands and give their rough classification. In section |IT1| we summarize our 
numerical results on the effects induced by resonant scattering in a simple fee crystal of 
dielectric spheres. Contrary to the suggestions made previously in the literature 0, no 
spectacular effects may be expected. Finally, in section |V| our conclusions are presented. 



II. RESONANCE SCATTERING AND BAND STRUCTURE OF PHOTONIC 

CRYSTALS 



Multiple scattering of classical waves in the presence of resonant scatterers was first 
studied in a disordered medium. There resonances cause a scattering delay due to the 
storage of wave energy inside a single scatterer, resulting in a sharp decrease (depending on 
the filling fraction) of the transport velocity ve of light |TB|. The discussion of resonance- 
induced effects in an ordered medium has been initiated by John 0. He suggested that 
the coinciding resonance and Bragg scattering is the most favourable condition for opening 
a gap in the spectrum. In the case of spheres this leads to / = 1/(2^/6) p[. Later on, 
Zhang and Satpathy Jl^ noticed that a pseudogap in the band structure of an fee lattice of 
dense spheres corresponds to a Mie resonance. Recently, Ohtaka and Tanabe 
KKR method described earlier in 



18 1, using the 
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made an attempt to relate the Mie resonances to the 
photonic bands. They paid attention to flat bands which correspond to "heavy photons" in 
analogy to their electronic analogue, heavy fermions. For heavy photons, the group velocity 
can be as low as ~ c/100 The role of Mie resonances in bonding of spheres in photonic 



crystals was investigated by Antonoyiannakis and Pendry 



In our paper, we shall attempt to investigate the complementary relation to that discussed 
by Ohtaka and Tanabe |18], namely, to what extent the single sphere resonance frequencies 
are related to band gaps, a question asked first by John [||. In obtaining his result, John 
0] used a simplified picture in which, regardless of the propagation and polarization, the 
photon always encounters precisely the same periodic e{r) variation, resulting in a one- 



dimensional Kronig-Penney's (KP) model |2y] . However, this approximation is only partially 
justified in the short-wavelength limit when light in a sphere behaves like a field in a slab 
with thickness 2rs. For real crystals the situation is different and the above coincidence 
condition can only be met in restricted regions on the surface of the Brillouin zone. Also, 
with increasing dielectric contrast 6, dispersion curves lv = ui(k) become nonlinear and 
the diffraction condition is modified as compared to the Bragg case and fulfilled at lower 



frequency [T^. Moreover, the model neglects polarization dependent effects. 

However, it is quite difficult to relax the simplifying assumptions of the one-dimensional 
KP model made by John 0. In our discussion of resonance-induced effects in photonic 
crystals we shall employ the machinery of the on-shell multiple-scattering theory (MST) 
(see for electrons and p2| for photons). Note that the on-shell MST is already required 
in the one-dimensional KP model if one wants to go beyond the dispersion relation and 
obtain Green's function or the local density of states p3|. The unique feature of the on- 
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shell MST is that, for nonoverlapping (muffin-tin) scatterers [pl| , p4| (the present situation), it 
disentangles single-scattering and multiple-scattering effects (see ||2^ for a recent discussion). 



For example, the total T-matrix per unit cell can be written in the form 

T = 1/{R-^ - B) = R/{1 - BR). (1) 

Here all quantities are ordinary matrices. R stands for the on-shell reaction matrix (also 
known as the k-matrix [^) of the scattering sphere (of all scatterers inside the lattice 
unit cell in the case of a complex lattice). It is diagonal in the angular- momentum basis 
and can be written as Rll' = —^ll' tan 17^/0", where rj^ is a phase shift. Here L is a 
composite index which labels all the spherical harmonics in the irreducible representation 
of the rotation group characterized by the principal angular-momentum number I and, in 



the case of electromagnetic waves, it carries an additional polarization dependent index |^ 
B = B{(T, k) in Eq. (|l|) is a matrix of so-called structure constants which accounts for the 
periodicity of the lattice. It depends on a and the Bloch momentum k. The R-matrix is 
singular at Mie resonance frequencies and i?((T, k) is singular whenever cr^ = (k + K„)^, 
where K„ is a vector of the dual lattice. 

The exact eigenmodes of a crystal are determined by poles of the total T-matrix, the 
latter being the zeros of the determinant of a hermitian matrix R~^ — B, 

det(i?"^ -B)=0. (2) 

Equation (^ is the familiar Korringa-Kohn-Rostocker (KKR) equation [^] in band structure 



theory. Recently, we have succesfully used this approach to calculate the band structure of 
electromagnetic waves in a simple fee lattice of dielectric spheres 0, to establish a simple 
analytic formula describing the width of the lowest lying stop gap in the (111) crystal 
direction (the L point of the Brillouin zone p8|) [^], and to calculate the properties of the 



local DOS in one dimension ||23|. In general, the higher frequency, the higher the value 
of ^max is needed to ensure convergence. In order to reproduce the first band and the 
linear part of the spectrum, Imax = 1 is enough. In general, the size of a secular equation 
is reduced by almost a factor 10 compared with that in the plane-wave method [p|,|7|,|8|JT7[| 
which customarily requires well above a thousand plane waves. The precision of the elements 



of the secular equation is determined by the standard Ewald summation which yields 
structure constants up to six digits. 

At first sight, there seem to be large differences in band formation between the electrons 
and electromagnetic waves. In the case of electrons, in the medium between two atoms, waves 
are evanescent in nature, while electromagnetic wave propagates unattenuated between two 
scatteres. However, electrons are strongly interacting with each other. If the interactions are 
taken into account, an effective-Hamiltonian single-electron picture emerges of a near-free 
electron with a positive energy which moves unattenuated between two scatteres like the 
electromagnetic wave does [^. Therefore, the same principles apply to the classification 
of the electronic and the photonic bands. For electrons, a rough classification of the bands 



can be obtained if the singularities of R and B are well separated EH]. Thus if -R(cr) is 



sufficiently small near a singularity of B{a, k), where a = yE and E is the electron energy, 
a band is formed with the dispersion relation 

o;(k) ^ |k + K„| +i?(a). (3) 
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In more intuitive terms, the formation of such a band results from the formation of standing 
waves in a crystal and it is appropriate to call such a band the Bragg band. The second type 
of band can form near a singularity of R at cuq, R = r/(co'^ — uJq). This is, for example, the 
case in the transition and noble metals. If F and B{a, k) in the vicinity of lvq are sufficiently 
small, a (usually very narrow) resonance band is formed with the dispersion relation 

u{k)^[u;o + TB{ao,k)]'^'- (4) 

Eq. is in agreement with the observation made in [|T^ that the resonance-band width 
is comparable with the lifetime broadening of the DOS profile for a single sphere. In more 
intuitive terms, formation of the resonance band can be understood as resulting from the 
broadening of individual resonances when they start to feel the presence of each other, similar 
to the formation of the electronic bands from individual atomic levels in the tight-binding 
limit |]2B|. Such usually very narrow resonance bands describe "heavy photons'' 



This shows that one can associate with a resonance a (resonance) band. Is it possible to 
associate with a resonance a gap? The answer is yes. However, a hybridization of bands must 
take place. We speak of hybridization if the singularities of R and B cannot be well separated 
and neither a pure Bragg nor a pure resonance band is formed. Under certain conditions 
the two bands can hybridize in such a way as to create a gap over the approximate energy 
range of the original unhybridized resonance band. An example of such a hybridization 
is provided by transition metals with characteristic d (or /) resonance which couples with 



extended band states by tunneling In the latter case a broad s — p band hybridizes 
with a narrow d band in such a way as to create a gap over the approximate energy range 
of the original unhybridized d resonance band ||30|| . 

In the following section, we investigate the hybridization of photonic bands in an fee 
lattice of dielectric spheres. Although the same principles apply to the classification of the 
photonic bands as to those for electrons, is does not mean that the respective band structures 
are qualitatively similar. There is significantly difficult to open a gap in the spectrum of 
electromagnetic waves than in the case of electrons. Moreover, a gap often does not open 
between the lowest lying bands, as in the case of electrons, but in an intermediate region. 
The origin of this difference lies in a different behaviour of R and B and can be rather easily 
understood. Indeed, in the tight-binding picture of band formation [^, individual bands 
results from the broadening of corresponding atomic levels when the atoms start to feel the 
presence of each other. The largest gap between atomic levels is between the lowest-lying 
energy levels. Therefore, for a lattice of atoms, one expects to find a gap essentially between 
the first and the second energy band, with the gap between higher bands scaling down to 
zero. On the other hand, for a dielectric scatterer and Maxwell's equations, bound states 
are absent. They are replaced by resonances. Moreover, if the wavelength is small compared 
to the size of the spheres, one can use geometric optics, while in the opposite limit of long 
wavelengths, the Rayleigh approximation applies. In neither case does a gap open in the 
spectrum. Therefore, if a gap is present in the spectrum, it should be in the intermediate 
region between the two hmiting cases (see, however, the case of a diamond lattice ( 0, figure 
2), which is a complex lattice). The same applies to the localization of light which is also 
expected at some intermediate frequencies. 

Since opening of a gap in the spectrum of electromagnetic waves in much difficult than 
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in the electronic case, one expects also that hybridization will be weaker in the former case. 
These expectations are confirmed in the following section. 

III. RESULTS FOR AN FCC LATTICE OF DIELECTRIC SPHERES 

In our case of an fee crystal of homogeneous spheres, we looked for a correspondence 
between Mie resonance frequencies and (i) the lowest lying stop gap in the (111) crystal 
direction (the L point of the Brillouin zone [^) in the case of both "dense" spheres (e^ > Eb) 
and "air" spheres (e^ < Sb), (ii) the full gap in the case of "air" spheres. (Note that there 
is no full gap in the case of "dense" spheres 0,^.) Apart from numerous experimental data 
now available [0,^, there are at least two other reasons to chose the L-gap. First, the width 
of the first stop gap often takes on its maximum at the L point and, second, experimental 
techniques make it possible to grow collodial crystals such that the L direction corresponds 
to normal incidence on the crystal surface. 

In the case of "dense" spheres we found that for sufficiently high dielectric contrast 6 the 
L-gap can be associated with the IMl resonance. (Our notation lAn for a Mie resonance 
is such that / is the angular momentum, A stands for either electric (E) or magnetic (M) 
mode, and n is the order of the resonance with increasing frequency in a given I A channel.) 
For all filling fractions one finds that as 6 increases over a critical value (see Figure [^), the 
IMl resonance "descends" from above to the L-gap and stays inside it close to the midgap 
frequency. We verified this behaviour for 6 up to 100. 

In the case of "air" spheres the hybridization of bands is much less pronounced compared 
to the case of "dense" spheres. We observe a correspondence between the L-gap and a Mie 
resonance only for particular filling fractions and dielectric contrasts. Filling fraction f ^ 0.5 
or higher is required for hybridization to occur. However, for a given /, the required dielectric 
contrast for the onset of hybridization can be as little as half the value for the dense sphere 
case. For / = 0.6 one finds the lEl resonance trapped inside the L-gap already for 6 > 8. 
In the close-packed case the lowest lying lEl resonance descends actually below the L-gap 
and, if (5 > 16, a hybridization occurs at the frequency corresponding to the second, IMl, 
resonance. 

Hybridization in the case of the lowest L-gap is only partial since the gap does not 
extend over the whole Brillouin zone. Let us, therefore, look at the full band gap which 
can be opened only in the air sphere case (one can open just a single full gap here 0J^). 
In contrast to that of the L-gap, the opening of the full gap requires a certain threshold 
dielectric contrast which rapidly increases as / decreases from the close-packed case [^,|^ 
(see figure 0). Hybridization follows the irregular pattern seen in the case of the L-gap. For 
example, for / = 0.6 hybridization occurs if 5 > 36 where the 2M1 resonance descends from 
above to the full gap, soon followed by the 3E1 resonance (2M1 < 3E1). Both resonances 
seem to be locked inside the full gap (at least up to 5 = 100). Interesting behaviour is 
found for / = 0.64. First, the full gap opens at 5 ~ 10.6 between the closely lying 2M1 and 
3E1 resonances. As Sb increases, the 3E1 resonance moves across the gap and stays closely 
below the lower edge of the full gap. Hybridization around the 3E1 resonance occurs only 
for 6 = (11.5, 13). Outside this interval of 6 no resonance is inside the full gap. The case in 
which / = 0.68 shows an anomalous behaviour: no resonance frequency is inside the full gap 
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or close to the the full gap edge. On the other hand, in the close-packed case (/ ~ 0.74) one 
can find three closely lying resonances, namely, 3M1, 4E1, and 1E2 (3M1 < 4E1 < 1E2), 
inside the full gap. 3M1 descends to the full gap around Sb = 12, soon followed by the 4E1 
and 1E2 resonances. For Sf, > 16 all three resonances are already inside the full gap and 
seem to remain trapped there (at least up to Sb = 100). 

Apart from the hybridization, another intriguing issue is that of whether the width of a 
gap can be enlarged due to the Mie resonances as the one- dimensional KP model suggests 
0. Once a resonance was inside a gap we did not find any significant effect on the gap 
width. A good illustration of this is the behaviour of the relative width A*" of the full gap 
(the full gap width divided by the midgap frequency) at / = 0.68 on the one hand and 
/ = 0.6, 0.64 and 0.74 (close-packed) on the other hand (see figure |). In the first case, the 
gap does not correspond to any Mie resonance, while in the second case, there are up to 
three Mie resonances inside the gap. Had there been some effect of a Mie resonance on the 
gap width, one would have observed either a sudden increase in A*" after hybridization sets 
in, or an anomalously small gap width for / = 0.64 (if 6 > 13) and / = 0.68. Instead, one 
sees a monotonical increase of A*" as 6 and / increase. However, something different may 
happen if a Mie resonance frequency is close to a gap edge, before hybridization sets in. 
Under certain circumstances one can observe a resonance-induced widening of a relative gap 
width by up to 5%. Figure |^ shows the relative L-gap width A^ for air spheres and / = 0.6. 
Around 6 = 7.99 the lEl resonance crosses the upper edge of the L-gap, which results in 
a local enhancement of A£ by ~ 5%. As a function of 6, this widening occurs within the 
very narrow interval of 5 G (7.986, 7.998) where it shows a flat peak. In the air-sphere case 
one can be sure that this widening of a gap can be entirely attributed to a resonance. As 



shown recently ||2^, the lowest lying L-gap for an fee lattice of dielectric spheres can be 
understood in terms of simple quantities, namely, the volume averaged dielectric constant, 
e = fss + (1 — f)eb, the volume averaged e^(r), = [fel + (1 — f)el], and the effective 



dielectric constant Seff- The effective dielectric constant can be well approximated PJ5^ by 



Maxwell-Garnett's formula 33 



e,ff^eb{l + 2fa)/{l-fa) (5) 
where, for a homogeneous sphere, the polarizability factor 

a= {ss- eb)/{es + '2^eb). (6) 
Then the absolute L-gap width is approximated within 3-6% (depending on /) by the 



formula 29 



AL=C(/)(v^-ee//)'Ye (7) 

where 

C(/)=Co + 0.14/(2/„-/)//^. (8) 

Here Cq ~ 0.74 is the minimal value of C and fm is the fllling fraction for which C{f) takes 
on its maximal value. C{f) takes on its minimal value Co at the extreme fllling fractions 
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/ = and / = 0.74, and its maximal value is Cm ~ 0.88 at fm ~ 0.74/2. The factor 0.14 in 
the interpolation formula (§) is the difference Cm — Cq. Let ki be the length of the Bloch 
vector at the L point of the Brillouin zone. In units where the length of the side of the 
conventional unit cell of the cubic lattice is A = 2, one has ki/n = -\/0.75. The formula 

= 2Tin,ff/\L/kL (9) 

then describes for 1 < 5 < 100 with the relative error ranging from 4% for / around 
0.2, to the relative error ^ 8% for the close-packed case. Formula (|^) describes in the 
air sphere case irrespective whether there is (/ > 0.5) or is not (/ < 0.5) a Mie resonance 
within the L-gap frequencies. However, the formula is violated in the narrow interval 5 G 
(7.986, 7.998) where the local widening of A^ takes place. 



IV. CONCLUSION 

Using the photonic analog of the KKR method P,p2|, we have investigated Mie- 
resonance-induced effects for an ordered medium, namely, for a simple fee lattice of homo- 
geneous spheres in three dimensions. Our work partially fills the gap in the understanding 



of resonance-induced effects compared to a disordered medium |]T6[. We showed that the 
same principles apply to the classification of both electronic and electromagnetic bands in 
a periodic medium, although their respective qualitative behaviour may be different. For 
example, due to the absence of bound states for a single scatterer, it is much more difficult 
to open a full gap in the spectrum of electromagnetic waves than for the case of electrons. 
Still, under certain conditions, one can identify a resonance band, and a hybridization of 
the Bragg and resonance bands can take place, leaving behind a gap over the approximate 
energy range of the original unhybridized resonance band ||30|. We investigated both a par- 
tial hybridization in the case of the first stop gap in the (111) crystal direction (the L-gap) 
and a full hybridization in the case of the full band gap. For dense spheres {Ss > e^), partial 
hybridization around the lowest Mie resonance occurs for all filling fractions /, once the 
dielectric contrast 6 reaches a critical value 5c(/) (see Figure |1|). For air spheres (e^ < Eh), 
the hybridization follows an irregular pattern and it can be observed only if / > 0.5. How- 
ever, the value (5c(/) can be as little as half the corresponding Sc{f) in the dense-sphere 
case. Near close packing (/ ^ 0.74), the partial hybridization occurs around the second Mie 
resonance. The full hybridization in the case of air spheres (there is no full gap in the dense 
sphere case [0-0|) follows again the irregular pattern seen in the case of partial hybridization. 
Sometimes the full gap opens in a frequency region which does not correspond to any Mie 
resonance (/ = 0.68), but in other situations there are up to three Mie resonances in the 
frequency range corresponding to the full gap (/ = 0.74). The resonant-induced widening of 
a gap can occur if a Mie resonance is about to cross the edge of the gap (see Figure |^). The 
values of the dielectric contrast required to observe some of these effects are within exper- 
imental reach at microwaves and, in the air-sphere case, even at optical and near-infrared 
frequencies. However, unlike in the case of the one dimensional Kronig- Penney model 0, 
we did not find any evidence that if some of Mie resonance frequencies fall inside a gap, 
then this leads to its significant widening. Contrary to the suggestions made previously in 
the literature pi, no spectacular effects may be expected. This is probably related to the 
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known fact that the hybridization of bands in higher dimensions is weaker pO[. Also, since 
the opening of a gap in the spectrum of electromagnetic waves in much difficult than in the 
electronic case, it is natural that also hybridization effects are weaker in the photonic case. 
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FIG. 1. The value of the dielectric contrast at which the lowest Mie resonance enters the L gap 
versus filling fraction in the case of an fee lattice of homogeneous dense spheres. 
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FIG. 2. The relative width of the full gap (the full gap width divided by the midgap 
frequency) in the case of an fee lattice of homogeneous air spheres as a function of the refractive 
index contrast \/5 for different filling fractions /. 
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FIG. 3. The relative L-gap width (the L-gap width divided by the midgap frequency) for 
an fee lattice of air spheres with the sphere filling fraction / = 0.6 plotted against the refractive 
index contrast \/6. Around 5 = 7.99, the lEl resonance crosses the upper edge of the L-gap which 
results in a local enhancement of the relative L-gap width, which is shown as a very narrow peak 
in the main figure. The inset shows the detailed view of the peak. 
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